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Evaporation of ices near massive stars: models based on laboratory 
TPD data 
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ABSTRACT 

Hot cores and their precursors contain an integrated record of the physics of the collapse 
process in the chemistry of the ices deposited during that collapse. In this paper, we present 
results from a new model of the chemistry near high mass stars in which the desorption of each 
species in the ice mixture is described as indicated by new experimental results obtained under 
conditions similar to those hot cores. Our models show that provided there is a monotonic 
increase in the temperature of the gas and dust surrounding the protostar, the changes in the 
chemical evolution of each species due to differential desorption are important. The species 
H2S, SO, SO2, OCS, H2CS, CS, NS, CH3OH, HCOOCH3, CH2CO, C2H5OH show a strong 
time dependence that may be a useful signature of time evolution in the warm-up phase as the 
star moves on to the Main Sequence. This preliminary study demonstrates the consequences 
of incorporating reliable TPD data into chemical models. 

Key words: stars: formation; ISM: abundances; ISM: clouds; ISM: molecules 



1 INTRODUCTION 

Conventionally, the first indirect manifestation of a young mas- 
sive star is provided by the presence of 'hot cores'. Hot cores are 
small, dense, relatively warm, optically thick, and transient ob- 
jects detected in the vicinity of newly formed massive stars. It 
is well-known that they exhibit a range of molecular species un- 
like those found in quiescent molecular clouds: in comparison, hot 
cores are richer in small saturated molecules and in large organic 
species (e.g. Walmsley and Schilke 1993). The characteristic hot 
core chemistry is believed to arise from the evaporation of the icy 
mantles that accumulated on dust grains during the collapse that 
led to the formation of the massive star; solid-state processing in 
the ice and gas-phase chemistry after evaporation also occur (cf. 
Millar 1993; Nomura and Millar 2004). Moreover, recent observa- 
tions do indicate the existence of possible precursors of hot cores 
(e.g Molinari et al. 2002; Beuther et al. 2002). Thus, hot cores and 
their precursors may be regarded as containing an integrated record 
of the physics of the collapse process in the chemistry of the ices 
deposited during that coUapse.It is the task of astrochemistry to un- 
ravel this record from the observations. 

Viti & Williams (1999) noted that the duration in which the 
grains are warmed from very low (~ 10 K) temperatures to the tem- 
perature at which H2O ice desorbs and above is determined by the 
time taken for a pre-stellar core to evolve towards the Main Se- 
quence. This time is believed from observational and theoretical 
studies to be fairly brief, ~10''-10^ y (e.g. Bernasconi & Maeder 
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1996, hereafter BM96), and in some contexts this 'turn-on' of the 
star may be regarded as instantaneous. However, in the context 
of hot cores, which have a duration similar to that of the turn- 
on time of the hot star, the relatively slow warming of the grains 
must be taken into account. Evidently, highly volatile species will 
be desorbed from the ices at early stages, while the more refrac- 
tory species will only be desorbed at later stages. The chemistry of 
the hot cores must reflect this temperature-driven evolution, and the 
evolutionary chemistry should provide, in principle, a direct mea- 
sure of the stellar turn-on time. Thus, hot cores contain not only 
a record of the collapse process but also of the ignition history of 
the star. The de-convolution of hot core observational data should 
constrain both physical processes. 

The Viti & Williams study, in the absence of reliable data, 
made a gross simplification in the description of desorption from 
mixed ices; their paper emphasised that appropriate temperature 
programmed desorption (TPD) experiments on mixed ices under 
near-interstellar conditions were essential. Such definitive measure- 
ments have now been made at the University of Nottingham Sur- 
face Astrochemistry Experiment. Eraser et al. (2001) reported a 
TPD study of pure H2O ice, and CoUings et al. (2003a) examined 
the TPD of H2O/CO ice mixtures. Both experiments provided data 
that is in striking contrast to the assumptions made in models of hot 
cores. In particular, CoUings et al. showed that CO desorbed from 
an H2O/CO ice mixture occurs in four temperature bands, rather 
than the assumed CO desorption at temperatures above about 25 K. 
They also found that the proportion of the CO desorbing at each 
phase depended on the total CO concentration and on the prepa- 
ration of the sample. These laboratory results suggest that in hot 
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cores the CO being injected into the gas phase from the ices should 
occur in four phases, rather than a single pulse. 

The chemistry of interstellar ices is of course much more com- 
plex than a simple H2O/CO ice mixture, although H2O is always 
the dominant partner (Whittet 2003). Therefore, similar TPD ex- 
periments on a variety of ices have been carried out at Notting- 
ham (Collings et al. 2004). In addition, the properties of some other 
species with respect to desorption have been considered. It has been 
found possible to categorize all species relevant to hot cores as 
H20-like (i.e. behaving as determined by Fraser et al. 2001); CO- 
like (as described by Collings et al. 2003a), or placed in several 
other intermediate categories. Thus, it is now possible to attempt 
to model the chemistry near high mass stars in which the desorp- 
tion of each species in the ice mixture is described as indicated by 
experimental results obtained under conditions similar to those hot 
cores. This paper reports the results of such a study. 

We model the chemistry of gas and the deposition of the ices 
during the collapse, and subsequently the desorption and chemistry 
in the gas phase after the star begins to warm the core. Our purpose 
is to demonstrate the consequences of incorporating reliable TPD 
data into chemical models. Specific sources will be modelled in 
later papers. 

In Section 2 we describe the new TPD results. Section 3 gives 
descriptions of the physical and chemical aspects of the model, and 
the results are presented in Section 4. Section 5 gives a discussion 
of the results and conclusions. For simplicity, throughout the paper, 
we shall use the word 'hot core' to include any pre-stellar core such 
as those observed by Molinari et al. (2002 and references therein) 
which are believed to be possible hot core precursors. 



2 EXPERIMENTS 

In the accompanying paper, Collings et al. (2004) report results of 
TPD experiments for a large number of molecules relevant to hot 
core chemistry. Details of these experiments are not discussed here. 
Molecular were assigned to one of three several categories of des- 
orption as follows: (I) CO-like; (2) H20-like; and (3) intermediate. 
The work also enables us to make inferences about likely behaviour 
of some species not so far studied, and in doing so we have created 
two further categories: (4) reactive; and (5) refractory. For each cat- 
egory we have used the information in Collings et al. (2004), and 
previous studies of the CO/H2O system (Collings et al. 2003a,b) to 
estimate the fraction of a particular molecular species that is des- 
orbed in the various temperature bands. These desorption bands are 
from (i) the pure species; (ii) a mono-molecular layer on H2O ice; 
(iii) the amorphous to crystalline H2O ice conversion (the 'volcano' 
effect); (iv) co-desorption when the H2O ice desorbs; and (v) des- 
orption from the bare grain surface. The fraction of the species des- 
orbing in each of the five desorption processes is shown in Table I, 
assuming a relatively thick ice layer of 0.3 yum, as used in the TPD 
experiments. The desorption processes are tabulated in sequence, 
and the quoted value expresses the fraction of the species remaining 
in the ice at that time that is desorbed in the process. In the layered 
model of ice mantles, the species CO, N2 and O2, are thought to be 
found predominantly in an outer layer non-hydrogenated (apolar) 
ice, separate from the water dominated inner layer of hydrogenated 
ice (polar). As the mantle warms from 10 K, the apolar layer is 
able to diffuse into the porous structure of the water ice. We are yet 
to accurately measure the relative concentration at which the water 
will become saturated, but we estimate that this occurs at roughly 
0.25. Since the total concentration of CO, N2 and O2 is our model is 



0.35-0.4 relative to water, desorption of these species from the non- 
hydrogenated outer layer will be relevant, and we list 35% of each 
as desorbed from the solid phase. Methane should only be found 
within the hydrogenated inner layer, and we assume that this may 
also be true of NO and CS. Based on the CO results of Collings et 
al. (2003b), about 70% of molecules of CO-like species in the hy- 
drogenated layer will undergo mono-molecular desorption in hot 
core conditions, and the remainder will become trapped. The TPD 
results (Collings et al. 2004) suggest that trapped molecules are re- 
leased in volcano and co-desorptions of equal sizes for the species 
initially present in the water ice layer, and with a volcano desorp- 
tion twice the size of the co-desorption for the species that have dif- 
fused into the water ice. The water-like species were found to be re- 
leased entirely in the co-desorption process. To the experimentally 
tested species we have added polar organic molecules and other hy- 
drophilic molecules. To the intermediate category we have added 
untested non-polar organic species and other small molecules that 
we anticipate will interact with have moderate interaction with wa- 
ter ice. All of the intermediate species are assumed to be found in 
the hydrogenated ice layer. The TPD results (Collings et al. 2004), 
suggest that desorption will be roughly evenly split between the 
volcano and co-desorption events. However, for the smaller species 
C2H2, and H2S, a small fraction (~ 10%) is able to escape and un- 
dergo mono-molecular desorption. We anticipate that HCN, HNC, 
NS, and HCS will also behave in this manner. The nine species we 
have categorised as reactive are all ions or free radicals in the gas 
phase. It is not clear how they will behave when adsorbed in the 
grain mantle at 10 K, however as the ice warms they would be ex- 
pected to adopt a form quite different from their gas phase structure 
due to strong interactions with the water ice and other molecules. 
As with the water-like species, we expect that they will only escape 
to the gas phase in a co-desorption process. The refractory species 
are those such as atomic magnesium and S2 which are likely to 
be stable in the mantle but are strongly adsorbed and non-volatile. 
We intuitively expect that they will remain in the mantle as other 
species, including water, desorb around them. Unless some explo- 
sive desorption mechanism operated, we anticipate that they will 
'sink' through the desorbing mantle and become bound to the dust 
grain surface. Desorption from the dust grain surface is therefore 
the only relevant process for this category. 



3 THE CHEMICAL AND PHYSICAL MODEL 

The model employed is similar to that used in VW99 and Viti et al. 
(2001). For this study, we performed time-dependent single-point 
calculations for Av of ~ 600 mags as derived from hot core obser- 
vations. For simplicity we have not used a multi-point code, as we 
did in VW99. Our model is a two phase calculation: the first phase 
starts from a fairly diffuse (~ 300 cm"^) medium in atomic form 
(apart from a fraction of hydrogen in H2), and undergoes a free- 
fall collapse until densities typical of hot cores are reached (~ 10^ 
cm"^). During this time, atoms and molecules from the gas freeze 
onto the grains and they hydrogenate where possible, as in VW99. 
Note that the advantage of this approach is that the ice composi- 
tion is not assumed but it is derived by a time dependent computa- 
tion of the chemical evolution of the gas/dust interaction process. 
Also the abundances that we derive for the high density core are 
directly computed from the initial cosmic abundances. The main 
differences here with the VW99 models are i) the initial abun- 
dances which we corrected to match the latest findings by Sofia 
et al (2001) and are 1.0, 0.075, 4.45xl0-\ \.19x\Q-\ 8.52x10^^ 
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Table 1. Assignment of each molecular species into a category of desorption. The first column indicates the category, column 2 lists the species belonging 
to each category while columns 3-7 indicates the fraction of each species (remaining after the previous desorption) that desorbs at each desorption event. We 
have used italics for the expeiimentally tested species. 



Species 



Solid 
at T~ 20K 



Mono 
atEi 



Volcano Co-des Surface 



co-like CO, Ni, O2 

CS, NO, CH4 
Intermediate HCN, HNC, H2S, NS, HCS,C2//2 

0CS,CHiCN, C3H4,C02, SOi, HC3N,C2//4 
H20-like H2O. CH3OH. SO, CH2CO, C2H50H,A'//3,N02, H2CO,HCOOCH, H2CS 

Reactive HS2, CoH, OCN, O2H, C2H5, H2CN, HNO, C2, C2H3 

Refractory Mg, S2 



0.35 


0.7 


0.667 


1.0 








0.7 


0.667 


1.0 








0.1 


0.5 


1.0 











0.5 


1.0 














1.0 














1.0 

















1.0 



1.43xl0-^ 5.12xl0-^ respectively, for H, He, O, C, N, S, Mg; ii) 
the use of the UMIST rate 99 database (www.rate99.co.uk) rather 
than UMIST96; iii) the low initial density (cf. 10* cm^^ in VW99). 
The unusually high initial density in VW99 before collapse was 
employed for ease of computation. Our present models start from 
a more realistic low density diffuse medium and this has conse- 
quences for the final percentage of gas in the mantles. In VW99, 
effectively 100% of the gas species was frozen onto the grains once 
the final density was reached. However, in Viti et al 2001, we ar- 
gued that it is possible that freeze out is never total; here we assume 
that 99% of the gas is frozen out, leaving 1% in the gas (a more 
realistic description since some gaseous CO is always observed, 
even in regions where no millimetre continuum is detected, see e.g. 
Molinari et al. 2000). This difference affects the very start of Phase 
2 of our computations. These differences do not affect at all our 
qualitative conclusions. We also performed a model where the ini- 
tial abundances and rate file used matched the VW99 for an initial 
comparison (see Section 4.1). 

Phase 2 follows the chemical evolution of the remnant core. 
We simulate the effect of the presence of an infrared source in the 
centre of the core or in its vicinity by subjecting the core to an 
increase in the gas and dust temperature. The temperature reaches 
its maximum (~ 300K) at different times depending on mass of the 
new bom star (see VW99). 

Phase 2 is rather different from that in the VW99 models be- 
cause of the different treatment for the increasing temperature and 
the evaporation from the grains.In VW99 we adopted a linear de- 
pendence of the temperature with time. For this study we have 
opted for a more realistic approach where we derived the temper- 
ature of the gas as a function of the luminosity (and therefore the 
age) of the protostar. We used the observational luminosity func- 
tion of Molinari et al. (2000); in their Figure 6 the total emitted 
luminosity is plotted against the core mass for mass accretion rates 
of 10"' and lO"** solar masses yr"'. Using this plot we correlated 
the luminosity, and hence the effective temperature, with the age of 
the accreting protostar and found that a power law fitted the data. 
Our major assumption was then that the temperature of the gas and 
dust surrounding the accreting protostar increases according to the 
same power law as the stellar temperature, and we fitted this so that 
the maximum temperature reached by the gas at the contraction 
time, tc, i.e the time after which hydrogen starts burning and the 
star reaches the ZAMS; Table|2|lists the contraction times and the 
volcano and co-desorption temperatures as a function of the mass 
of the stars. The contraction times are taken from BM96. Note, that 
the precise contraction time and the exact behaviour of the tem- 
perature increase are mostly irrelevant to the purpose of this study: 
the only necessary assumptions are that i) the contraction time of 



Table 2. Contraction times (f^) and the volcano and co-desorption tempera- 
tures as a function of the mass of the stars. 



M 


tc 


Volcano Desorption 


co-desorption 


Mo 


10'' years 


Kelvin 


Kelvin 


60 


0.0282 


92.2 


103.4 


25 


0.0708 


90.4 


101.6 


15 


0.117 


89.5 


100.8 


10 


0.288 


88.2 


99.4 


5 


1.15 


86.3 


97.5 



a massive protostar is comparable to chemical timescales - all the 
evolutionary studies of the type performed by BM96 have shown 
that even the most massive stars (~ 120 Mq) contract in no less 
than 10,000 years and ii) the temperature increases monotonically 
and not following a step function. It is true that the volcano and the 
co-desorption temperatures decrease as a function of the slowing 
down of the heating rate so that the lower the mass of the star, the 
lower the volcano and the co-desorption temperatures (see Table|2j 
and Collings et al. 2004) however this decrease is very small and 
insignificant with respect of the evaporation times. 

Our new evaporation treatment is pseudo-time dependent in 
that we allow evaporation of a fraction of mantle species X (in a 
single step) when the temperature for a particular desorption event 
is reached. The temperature and fractions for each species are listed 
in Table 1. Note that the fractions are always expressed as fractions 
of the current amount and not the original. 



4 RESULTS 

The first task of this theoretical study is to determine whether the in- 
clusion of the five distinct desorption events occurring during TPD 
of dirty ice affect significantly the abundances predicted by previ- 
ous models, in particular the VW99 models. To isolate this effect 
we ran several models using the same initial conditions of VW99 i.e 
the same initial abundances (see Table 2 from VW99) and chemical 
database. We found that substantial differences are seen between 
these models and VW99. These differences are solely due to the 
new evaporation treatment. 

We then constructed a new grid of models using the revised 
initial elemental abundances listed in Section 3, and the latest 
database UMIST99. We ran several models covering stars with 
masses from 5 to 60 Mq. We have taken the contraction times 
from BM96 as the epochs at which the temperature of the core 
reaches 300K. The desorption temperatures were calculated by run- 
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Figure 1. Fractional abundances of selected sulfur-bearing species as a 
function of time for a 25 Mq star, equivalent to a contraction time of 70,000 
years. 



Figure 3. Fractional abundances of selected large molecules as a function 
of time for a 15 M© star, equivalent to a contraction time of 70,000years. 
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Figure 2. As in Figure 1 but for a 5 M© star, equivalent to a contraction 
time of 1.15x10'' yrs. 



ning simulations of water ice desorption (CoUings et al. 2003b) us- 
ing the power law temperature profiles fitted to the stars of each 
mass. Figures 1-3 show the fractional abundances obtained from 
our computations for selected species for respectively a M = 25 
Mg, M = 5 Mq, and M = 15 Mq star. 

The species that are most affected by the inclusion of the 
new evaporation treatment are sulfur-bearing species and the larger 
species. The sensitivity of sulfur-bearing species to physical and 
chemical variations during the lifetime of a hot cores is not surpris- 
ing as previous studies have shown that sulfur-bearing species are 
strongly time-dependent (e.g. Hatchell et al. 1998; Viti et al. 2001, 
Hatchell & Viti 2002). While larger molecules, such as CH3OH 



and CH2CO, maintain the same trends from a 5 to a 25 M© star, the 
behaviour of sulfur-bearing species depend on the heating rate and 
hence on the mass of the star. In particular, we find that the evap- 
oration of H2S at different evolutionary stages is one of the main 
key factor to the changes occurring in the sulfur bearing chemistry: 
for a 25 M© star the abundance of H2S itself varies by about one 
order of magnitude with time; as the mass of the star decreases this 
difference becomes larger, up to two orders of magnitude. H2S has 
two main peaks in abundance due to desorption from the volcano 
and co-desorption events, but each peak is followed by a quick and 
sharp decline due to dissociation. This fast and changing behaviour 
has several consequences on other sulfur bearing species; it is inter- 
esting to look at these changes in the ratios of selected sulfur bear- 
ing species e.g: for a 25 M© star, (i) H2S/SO2 is > 1 up to 50,000 
years, and < 1 for older cores, (ii) H2S/CS, H2S/H2CS, H2S/OCS > 
1 up to 150,000 yrs and < 1 for older cores; (iii) SO/CS, SO/H2CS, 
SO/OCS ~ 1 up to 65,000 yrs and < 1 after that. For a 5 M© star, we 
find similar trends but much more accentuated; for example, H2S, 
SO and SO2 show similar abundances for longer periods during the 
earlier phases of star formation, while H2S/SO2 is much less than 1 
after 190,000 yrs. In general H2S decreases much more drastically 
with time as we move toward the lower mass regime; this is be- 
cause, although the temperatures at which the evaporations occur 
are quite similar between the 5 and the 25 M© cases, the tempera- 
ture increases occurs more slowly for a 5 M© star, allowing more 
time for the gas H2S to be destroyed, before another H2S evapora- 
tion 'peak' occurs. The main destruction route for H2S is its reac- 
tion with HCO^; the latter becomes abundant once the temperature 
reaches ~ 21 K and the first CO evaporation peak occurs. From 
these trends, it is clear that the ratio of these sulfur-bearing species 
can be used to investigate the evolution of the early phases of high 
mass star formation, and that sulfur bearing species are therefore 
good 'chemical clocks'. 

Another set of species that can be used as evolutionary indi- 
cators are organic molecules. Figure 3 shows the time evolution for 
selected large species for a 15 M© star. From the figure, we note 
that large species, such as C2H5OH, HCOOCH3 and CH3OH, are 
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good indicators of old cores as tiiese strongly bound species are 
abundant in the gas only at late times. 

Although observationally less significant, we would also like 
to draw attention to other species (not shown) namely: CO varies 
by about 1 order of magnitude during the increase of the tempera- 
ture and therefore during the formation of a high mass star but, CO 
being a ubiquitous species, this difference, even if detectable, may 
not, by itself, be indicative. Finally, we find that N2H^ should be 
present in detectable quantities in young hot cores but disappears 
probably within the first 15,000 years of its lifetime, when the tem- 
perature has reached ~ 40K; this corresponds roughly to the mono 
evaporation of CO (~ 40 K) and of CH4 (~ 45 K). In fact the main 
destruction routes for N2H^ are its reactions with methane and CO. 
On the contrary ammonia should only be detectable in cores older 
than 30,000 years (for a 25 Mq) Intermediate aged high mass star 
formation cores should be in general nitrogen-deficient. 



5 CONCLUSIONS 

We have extended existing chemical models of high mass star for- 
mation cores by including the latest experimental results on des- 
orption from grains in the evaporation of icy mantles formed in 
star-forming regions. We find that distinct chemical events occur 
at specific grain temperatures. This supports our assertion that in a 
comprehensive model it should be possible to use these chemical 
events as a chemical clock tracing the ignition of the nearby star. 

We have shown that, provided there is a monotonic increase in 
the temperature of the gas and dust surrounding the protostar, the 
changes in the chemical evolution of each species due to differential 
desorption are important. Based on the recent TPD experimental 
results, the species that are most affected are: H2S, SO, SO2, OCS, 
H2CS, CS, NS, CH3OH, HCOOCH3, CH2CO, C2H5OH, and to a 
lesser extent CO, N2H+, NH3 and HCN, HNC. Our results imply 
that if a large enough sample of high mass precursors objects are 
observed in these molecules, we would gather enough information 
to determine their age and therefore the evolutionary stage of the 
high mass star formation core. 

Finally, we note that I) effects induced by the relatively slow 
warming of gas near stars are more pronounced near stars of mod- 
erate mass, since the turn-on time for stars of intermediate mass 
is much longer than for massive stars. Observations of cores simi- 
lar to hot cores are now being found for intermediate-to-low mass 
stars (e.g Cazaux et al. 2003); 2) the evolutionary trends found in 
this study depend on the temperature gradient. Hence, at the estab- 
lished 'hot core' stage, similar trends should be found as a function 
of distance from the central star, again due to a temperature gradi- 
ent (e.g Nomura and Millar 2004). 
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